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ABSTRACT 

We used the deep colour-magnitude diagrams (CMDs) of five rich LMC clusters (NGC 1805, NGC 1818, NGC 1831, 
NGC 1868, and Hodge 14) observed with HST/WFPC2 to derive their present day mass function (PDMF) and its 
variation with position within the cluster. The PDMF was parameterized as a power law in the available main-sequence 
mass range of each cluster, typically 0.9 ^ iti/Mq ^ 2.5; its slope was determined at different positions spanning from 
the very centre out to several core radii. The CMDs in the central regions of the clusters were carefully studied earlier, 
resulting in accurate age, metallicity, distance modulus, and reddening values. The slope a (where Salpeter is 2.35) was 
determined in annuli by following two distinct methods: 1) a power law fit to the PDMF obtained from the systemic 
luminosity function (LF); 2) a statistical comparison between observed and model CMDs. In the second case, a is a 
free input parameter in the CMD modelling process where we incorporate photometric errors and the effect of binarity 
as a fraction of unresolved binaries (/bin = 100 %) with random pairing of masses from the same PDMF. In all clusters, 
significant mass segregation is found from the positional dependence of the PDMF slope: a ^ 1.8 for R < 1.0 Rcorc and 
a ^ Salpeter inside i? = 2 ~ 3 Rcoro (except for Hodge 14, where a ~ Salpeter for 7? ~ 4 Rcore). The results are robust 
in the sense that they hold true for both methods used. The CMD method reveals that unresolved binaries fiatten the 
PDMF obtained form the systemic LF, but this efi'ect is smaller than the uncertainties in the a determination. For 
each cluster we estimated dynamical ages inside the core and for the entire system. In both cases we found a trend in 
the sense that older clusters have fiatter PDMF, consistent with a dynamical mass segregation and stellar evaporation. 
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1. Introduction 



The phenomenon of mass segregation in a stellar system, 
which means a preferential concentration of high-mass stars 
towards the centre and a preferential allocation of lower- 
mass stars towards the periphery, seems to occur in sys- 
tems with widely distinct physical properties. It has been 
detected and extensively studied in globular clusters (de 
Marchi & Paresce 1996; Andreuzzi et al. 2000; Howell et al. 
2000), open clusters (Raboud & Mermilliod 1998; Durgapal 
& Pandey 2001; Bonatto & Bica 2003, 2005), and even in 
star forming regions (Hillenbrand & Hartmann 1998; Stolte 
et al. 2002). 

A key role in the investigation of mass segregation was 
played by the Hubble Space Telescope (HST), which for the 
first time resolved the stars in the very centre of the rich 
star clusters in the Magellanic Clouds (MCs) . This provided 
a new "laboratory" for obtaining constraints on the physi- 
cal processes involved in star formation, where the possible 
universality of the initial mass function (IMF) is a central 
issue (Kroupa 2002). The main reason is that, unlike the 
Galaxy, the MCs present a great variety of clusters, includ- 
ing young and rich star clusters. 

There are two possible and distinct origins for this ef- 
fect of mass segregation: dynamical and primordial. The 
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first one is caused by the dynamical evolution of the clus- 
ter, where the stars tend to reach the equipartition of ki- 
netic energy due to stellar encounters (Spitzer 1987; Binney 
& Tremaine 1987). Therefore, the high-mass stars decrease 
their velocities, sinking towards the cluster centre, while 
the low-mass stars speed up and take higher orbits on av- 
erage. In a simplified discussion, the characteristic time- 
scale of dynamical mass segregation in a stellar system is 
given by r?T.iow/'7ihigh times the two-body relaxation time 
{tri) (Spitzer 1987), where toiow and TOhigh are the lowest 
and highest masses in the cluster, respectively. This rela- 
tion indicates that the time when mass segregation occurs 
scales with tri and can be very short if rrihigh >> mioyj. An 
example of a very young ( ^ 2 Myr) system that presents 
mass segregation, which can be interpreted as a dynami- 
cal effect, is the Orion Nebulae Cluster (ONC) (Kroupa, 
Aarseth & Hurley 2001). 

It is also important to note that dynamical mass segre- 
gation combined with stellar evaporation may lead to the 
preferential loss of lower-mass stars. Since these stars are 
more likely populating the outermost regions of a cluster, 
they are more prone to be unbound due to their lower bind- 
ing energy. Therefore one could expect that stellar clusters 
would have flatter global PDMF as they become dynam- 
ically older. In fact this effect is observed in open clus- 
ters (Bonatto & Bica 2005) and Galactic globular clusters 
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(Baumgardt & Makino 2003) , which are also modelled with 
N-body simulations. 

On the other hand, primordial mass segregation may be 
a natural outcome of star-formation theory, since a proto- 
cluster with higher central density should have a greater 
probability of forming proportionally more high-mass stars 
in its centre. Some scenarios propose mass segregation at 
the onset of star formation through interactions among the 
protostars, since the collision probabilities increase with 
density (Bonnel & Davies 1998); alternatively, accretion 
rates are enhanced with the mass of the accreting proto- 
star (Behrend & Meader 2001). 

Regardless of the physical mechanism used to account 
for the origin of mass segregation, there are several tech- 
niques for diagnosing and quantifying the effect, the main 
ones being based on stellar statistics. By counting stars in 
different annuli, one may search for variations in radial pro- 
files as a function of stellar-mass range or for changes in the 
slope of the luminosity function (LF) or of the present day 
mass function (PDMF) . Although there are uncertainties in 
the conversion of stellar luminosity into mass (de Grijs et al. 
2002a) , the second option provides a more direct constraint 
on N-body simulations that intend to recover the initial 
conditions of the cluster (Kroupa, Aarseth & Hurley 2001; 
Baumgardt & Makino 2003; Moraux, Kroupa & Bouvier 
2004 ) . Since the mass-luminosity relation is dependent on 
metallicity, its uncertainty may be efficiently reduced by 
precisely determining the cluster physical parameters, by 
means of careful modelling of its colour-magnitude diagram 
(CMD). A detailed enough CMD modelling should result in 
physical parameters predominantly limited by uncertainties 
associated to the models of stellar evolution, rather than to 
the data. 

The main goal of this work is to determine the spatial 
dependence of the PDMF slope of five rich Large Magellanic 
Cloud (LMC) clusters, namely NGC 1805, NGC 1818, 
NGC 1831, NGC 1868, and Hodge 14. In Kerber & Santiago 
(2005), we presented the analysis of deep CMDs from these 
clusters obtained with HST/WFPC2 in the F555W (- V) 
and F8I4 W (^ /) filters. Efficient use of the data was made 
by means of direct comparisons of the observed CMD (sta- 
tistically corrected for incompleteness and field star con- 
tamination) to model ones. By modelling the CMDs in the 
central region of each cluster we inferred the metallicity 
(Z), the intrinsic distance modulus ((to — M)q) and the 
reddening value {E{B — V)). We also determined the age 
(r) for NGC 1831, NGC 1868, and Hodge 14 

Santiago et al. (2001) and de Grijs et al. (2002ab), using 
the same data, analyse the mass segregation in these clus- 
ters by means of the spatial dependence of the LF slope. 
Their diagnostic was clear: all clusters present mass segre- 
gation, even the youngest ones (NGC 1805 and NGC 1818). 
Meanwhile, de Grijs et al. (2002ab) reach the same result 
by comparing the radial profile dependence with the stellar 
mass range. They also derive the PDMF and its variation 
with position within the cluster, but only for the youngest 
ones. Here we present this type of analysis for the five 
clusters, using the cluster parameters derived by Kerber 
& Santiago (2005). The PDMF slope (a) in different an- 
nuli was determined by following two distinct methods: 1) 
power law fit to the PDMF obtained from the systemic 
LF (hereafter LF method); and 2) a statistical compari- 
son (similar to Kerber & Santiago 2005) between observed 
and model CMDs (hereafter CMD method), where a is a 



Table 1. Main parameters of the clusters in the sample. 



Cluster A'cius Ri 



I Rn 



clus -* '•core -* "-max 

njn 



Z log(r/yr) (m - A/)o E(B - V) 



NGC 1805 2564 


5.5 


55 


0.008 


7.80* 


18.55 


0.03 


NGC 1818 3929 


10.1 


60 


0.004 


7.80* 


18.45 


0.00 


NGC 1831 7136 


18.3 


90 


0.012 


8.70 


18.70 


0.00 


NGC 1868 5675 


6.7 


70 


0.008 


8.95 


18.70 


0.00 


Hodge 14 1196 


7.4 


40 


0.006 


9.25 


18.55 


0.03 



*youngest isochrone available by Girardi et al. (2000). 



free input parameter in the CMD modelling process. The 
main difference between the two methods is that the CMD 
method potentially uses all the information contained in 
the CMD, including the effects of unresolved binaries and 
photometric uncertainties. 

The paper is divided as follows. In Sect. 2 we present a 
brief description of our data and the physical properties of 
the clusters. In Sect. 3 we present the two methods of de- 
termining the PDMF slope and their results. These results 
are then discussed in Sect. 4, where we also compare them 
with those available in the literature. Finally, in Sect. 5 we 
present a summary and our concluding remarks. 



2. The data 

We used data taken with HST/WFPC2 as part of a cycle 7 
project entitled "Formation and evolution of rich LMC clus- 
ters" (Beaulieu et al. 1999). For each cluster and a nearby 
field, images were obtained using the F555W (^ V) and 
F8I4 W (~ /) broad band filters. A detailed description of 
the photometry and sample completeness corrections can 
be found in Santiago et al. (2001) and Castro et al. (2001). 
Kerber & Santiago (2005) made a detailed study of the 
resulting CMDs in the central cluster regions, in order to 
infer the cluster's global parameters, such as age, metal- 
licity, foreground extinction, and distance. Their work was 
based on a detailed CMD modelling process. In brief, the 
modelling is based on the generation of synthetic CMDs to 
be compared to the observed one. For a realistic comparison 
to the models, the observed CMDs had to be corrected for 
selection effects, such as photometric incompleteness and 
field-star contamination. Random photometric uncertain- 
ties in the real data were also quantified and applied to the 
model CMDs. Finally, different statistical tools were used 
in the model vs. data CMD comparison. 

Table n lists the main parameters of the clusters in the 
sample. The cluster name is followed by the total number 
of stars (A^'cius) in its final CMD after being corrected for 
the effects mentioned earlier (see Kerber & Santiago 2005 
for details). This table also lists the cluster's adopted core 
radius (-Rcoro, as determined by Mackey & Gilmore 2003); 
and the radius at which the star density falls to the field 
density (i?max)- The metallicity (Z), logarithmic age (r), 
intrinsic distance modulus ((tti — M)o), and reddening value 
{E{B — V)) were kept fixed in the PDMF analysis and are 
consistent with the ones determined by Kerber & Santiago 
(2005). 
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3. Analysis 

The PDMF is considered here to be a power law: 



am) = 



dN 
dm 



Co m' 



(1) 



where ^o is a normahzation constant and a the PDMF 
slope. Different paramcterizations are found in the liter- 
ature, often making use of multiple power laws. Given the 
limited mass range of main-sequence stars in the CMDs 
studied here, use of a single slope is justified. 

To analyse the positional dependence of a within a clus- 
ter, we divided its CMD stars into subsamples, according to 
the distance from the cluster centre, defining several annuli. 
The inner and outer radial limits of each ring were chosen 
to ensure that each subsample would typically contain 600 
stars, allowing a statistically significant number of stars to 
contribute to each PDMF determination. 

To determine the PDMF slope (a) in each ring we used 
two distinct approaches, each one made up of several steps: 

1) derivation of the systemic LF; conversion of the LF into 
a PDMF, according to a given mass-luminosity relation; a 
power law fit to the PDMF obtained from the systemic LF 
We refer to this method as the LF method. 

2) generation of synthetic CMDs with different input 
PDMFs, keeping the binary fraction and model isochrone 
parameters fixed; statistical comparison between observed 
and model CMDs; determination of the best-fitting CMD 
models for each data CMD. In this case, the PDMF slope 
a is a free input parameter of the CMD models. As the 
statistical comparison in this case makes use of the infor- 
mation available in the entire CMD plane, we call it the 
CMD method. 

Both approaches require conversion from mass to lu- 
minosity (CMD method) or vice- versa (LF method). The 
mass-luminosity relation used in these conversions is pro- 
vided by Padova isochrones (Girardi et al. 2000) that are 
shifted in magnitude and colour by (m — M)o and E{B — V) 
and where the parameters (including Z and r) are as given 
in Table ^ When necessary, we interpolated the original 
Padova isochrone grid in metallicity (see Kerber & Santiago 
2005, Sect. 3.2 and their Fig. 13) in order to generate an 
isochrone with the quoted values. 

Each approach displays its own advantages and disad- 
vantages. The LF method, by construction, does not use 
all the two-dimensional information contained in the CMD 
plane. This may actually be considered an advantage, since 
the PDMF can be reliably recovered with a smaller number 
of stars, as each magnitude bin concentrates the informa- 
tion spread along the MS colour width. On the other hand, 
the effect caused by unresolved binaries may be crucial, 
since the CMD position of primary stars will be spread 
redwards and brightwards due to the presence of the secon- 
daries. If unaccounted for, unresolved binarity will cause 
some shallowing in the recovered PDMF, as the system 
masses resulting from their combined luminosities will be 
larger than the masses of the individual components. 

Several previous works have taken the effect of unre- 
solved binaries into account in the recovered PDMF or IMF 
from a systemic LF. Kroupa, Tout, & Gilmore (1991) have 
managed to reconcile the LF of local volume-limited sam- 
ples of Galactic field stars with the LF inferred from pho- 
tometric surveys of more distant by correcting them for 
the systematics caused by unresolved pairs. More recently. 



Kroupa (2001) showed that the single-star IMFs can be sys- 
tematically steeper by 0.5 between 0.1 < m < 1.0 Mq than 
the Galactic-field IMF. Studying LMC star clusters, Sagar 
& Richtler (1991) determined that the recovered PDMF 
from the systemic LF can become significantly steeper (by 
Aa ~ 1.0) if the binary fraction is large (/bin ~ 0.50) 
and a ^ 0.5. Similar techniques were applied by Sandhu, 
Pandey, & Sagar (2003) to correct the PDMF for interme- 
diate/old open clusters. 

In order to simulate the effect of unresolved binaries in 
PDMF determination we introduced it in the CMD mod- 
elling process. The well-known signature of unresolved pairs 
in a CMD (see Hurley & Tout 1998 for a demonstration 
of the effect) can be modelled in a straightforward man- 
ner by applying the combined fiuxes and colours from two 
stars to a given fraction (/bin) of systems. One caveat may 
be the uncertainties in the distribution (dN/dq) of sec- 
ondary/primary mass ratios (g = ?7^2/?7^l < 1.0, where 
mi and m2 are, respectively, the primary and secondary 
masses) . This may at first sight be regarded as an extra de- 
gree of freedom in the modelling process, since secondary 
star masses may not necessarily be drawn from the same 
distribution as the primary stars. However, this possibility 
poses a question of what a mass function is meant to be, as 
it would not be uniquely defined even in a single popula- 
tion. Here we adopt the assumption that secondary stars in 
binary systems have masses drawn from the same PDMF as 
primary stars or as single stars. To maximize the possible 
effect of unresolved binaries and therefore explore the most 
of this effect, we used /bin=100%. 

Furthermore, the CMD modelling naturally incorpo- 
rates the photometric uncertainties into the PDMF deter- 
mination, and can potentially incorporate other observa- 
tional effects that may influence the conversion of magni- 
tude and colour information into mass. The major disad- 
vantage of a CMD method is the prior PDMF parameteriza- 
tion, here modelled as a power law with one free-parameter 
for m > 0.80 Mq ( ^ observed lower-mass limit). Even 
though it does not directly affect the PDMF in the observed 
regime, we fixed a shallow slope a — 1.30 in the lower mass 
range (0.08 < tti/Mq < 0.80). This is consistent with the 
IMF proposed by Kroupa (2002) in this mass range. This 
relatively shallow slope in the low-mass range also yields 
an enhanced effect of unresolved binaries. Although in our 
simulations all stars have one companion, only pairs with 
q ^ 0.60 will significantly change the CMD position of the 
primary star. The fraction of such effective binary systems, 
/bin.cff, is '^ 20%. They are practically the only ones re- 
sponsible for the effect that unresolved binaries may cause 
in the CMD or in the LF. We refer to Tout (1991) for a 
demonstration of the sensitivity of PDMF slope with the 
mass-ratio distribution in the low-mass regime. 



3.1. LF method 

Mass segregation was observed by Santiago et al. (2001) in 
their study of the LFs of the clusters in our sample. Then 
de Grijs et al. (2002a) converted LFs into PDMFs and pre- 
sented their positional dependence for the younger clusters, 
NGC 1805 and NGC 1818. But those authors did not use 
CMD modelling techniques to investigate the behaviour of 
a as a function of distance from the cluster centre. Also, 
the mass-luminosity relations used in previous works were 
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Fig. 1. PDMFs at different concentric regions around NGC 
1805. The lines are power law fits to the PDMFs. The slope 
values and their corresponding fit errors are given in each 
case, as well as the inner and outer radii of each region. 



Fig. 3. Same as in Fig. [2 but now for NGC 1831 
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Fig. 2. Same as in Fig.[ll but now for NGC 1818 



not the same as those used here, which are based on care- 
ful constraints on the cluster's ages and metallicities, which 
come from our previous CMD modelling. 

Notice that the mass-luminosity relation is a key ingre- 
dient in conversion from an LF into a PDMF, a traditional 
procedure widely used in stellar statistics. By definition. 




0.1 0.2 0.3 0.1 0.2 0.3 

log(m/ma) logCm/mJ 

Fig. 4. Same as in Fig. H but now for NGC 1868 

the LF ($(A^555)) in the M555 absolute magnitude is given 

by 



$(M555) 



diV 
dM. 



555 



(2) 



where dA^ is the number of individual stars that have ab- 
solute magnitude inside (M555, M555 -I- dMsss). 

If the PDMF is expressed in linear mass bins (.^(m)), it 
is related to the LF as follows 



diV 
C(m) = — = -<i>(M555) 
am 



dm(M555) 
dM555 



(3) 
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Fig. 5. Same as in Fig. ^ but now for Hodge 14 



where the 771(^/555) is the mass-luminosity relation. Here 
one clearly sees the importance of the choice of suitable 
stellar evolution models, as they provide not only the mass- 
luminosity relation, but also its derivative. As previously 
discussed, it is important to keep in mind that the observed 
photometric data suffer from unresolved binarity and there- 
fore allow us to construct only the systemic LF. As a conse- 
quence, there is not a unique mapping from one M555 mag- 
nitude to one stellar mass; the LF method may be seen as 
a first attempt to recover the PDMF. We refer the reader 
to Sect. 14.11 where we present an approach to evaluating 
the effect of unresolved pairs in the recovered PDMF from 
a systemic LF. 

Analogous to the LF, the PDMF may also be expressed 
in logarithmic (logio) mass bins {^-^{m)), in which case it 
is related to £,{m) by 



Um) 



dN 
dlogm 



-^M 



dlogTTT, 



dm 



— loge m £,{in). (4) 



Also parameterizing the logarithmic PDMF as a power 
law, we have £,i,{m) = ^l,o "t- : whose slope F is correlated 
with a by r = —(a — 1). 

Taking the logarithm on both sides of Eq. Q] we have 
logf = -alogm + log^o or log^L = -(a - l)logm + log^o - 
0.362. In other words, in a log-log plot, both the linear and 
logarithmic power law PDMFs should behave as a straight 
line, whose slope yields a directly. 

Figures I1I5I present the behaviour of the PDMFs re- 
covered directly (without any treatment of the unresolved 
binarity) from the systemic LFs at different concentric re- 
gions around each cluster. For convenience, we determine 
^L and plot it in a logarithmic scale, where we shift the 
values of log(^L) up or down to avoid overlapping in the 
panels. In each case a was obtained by means of a linear 
fit, shown as a straight line in the figure panels. Panel a) 
shows the PDMF for all stars in each cluster and the PDMF 



for the stars within the core. In all cases a is smaller in the 
core than for the entire cluster, which indicates mass segre- 
gation. Besides, these figures show the expected signature 
of mass segregation: the PDMF is shallower (smaller a) in 
the inner regions (panel b) than in the outer regions (panel 
d). Also, in general, the single slope fit is a very adequate 
description of the data. The exceptions to this rule often 
result from a sudden drop at the low-mass end (log(m/M0) 
~ 0.0, m ~ 1.0 Mq) of the PDMFs in the central regions 
of the richest clusters, such as NGC 1868. These deviations 
from a power law are very likely due to residuals in the 
correcting for incompleteness effects. 

Figure IHI plots a as a function of distance R from the 
cluster centre. This distance is expressed in units of the core 
radius, listed in Tabled The effect of mass segregation is 
again clearly seen for all clusters. The vertical bars on each 
point represent the linear-fit uncertainties in a. The hori- 
zontal bars just show the limits of the annuli. Some con- 
centric regions are partially overlapping in order to yield 
a more continuous behaviour of the PDMF slope and to 
increase the statistical significance in each PDMF determi- 
nation. The dotted-horizontal lines constrain the la range 
of a fitted to the global PDMF. This range can be com- 
pared with the Salpeter (1955) value {a = 2.35), marked as 
a dashed-horizontal line. In all the others clusters, besides 
Hodge 14, a{R) has a Salpeter value around i? ~ 2 — 3 
Rcoro- In the outer regions, the relation a vs. i?(i?corc) flat- 
tens for NGC 1818, NGC 1831, and NGC 1868. This may 
be the result of the dynamical loss of lower-mass stars from 
the clusters. 

3.2. CMD method 

Only the systemic LF, which suffers from the effect un- 
resolved binaries, is directly extracted from the data, as 
discussed in the introduction to this section. Therefore, we 
deal in this section with the 2D CMD modelling process, 
an approach that is capable of explicitly taking unresolved 
binarity into account. 

Our CMD modelling process and the statistical tech- 
niques of CMD comparisons were extensively explained in 
Kerber & Santiago (2005) and Kerber et al. (2002). Here 
we only underline some important aspects. The modelling 
process assumes that the cluster is a single stellar popula- 
tion (SSP) that generates a synthetic main-sequence (MS) 
in the CMD plane, where we introduce as model inputs 
the information about metallicity (Z), age (r) (given by 
a Padova isochrone; Girardi et al. 2000), intrinsic distance 
modulus ((m - M)o), reddening value {E{B - V)), PDMF 
slope (a), and fraction of unresolved binaries (/bin)- By ex- 
ploring a regular model grid, we may then find the best 
models by means of statistical comparisons carried out in 
1 and 2 dimensions and, therefore, the physical parame- 
ters that best constrain the cluster CMD. By modelling 
the CMDs in the central region, Kerber & Santiago (2005) 
inferred Z, (m — M)q, and E{B — V) for each cluster. They 
also determined the age of three clusters: NGC 1831, NGC 
1868, and Hodge 14. A set of values consistent with these 
determination is listed in Tabled 

We here model the CMDs separated by annuli in order 
to investigate the positional dependence of the only remain- 
ing free parameter: a. Figure[7|illustrates the CMD method 
for NGC 1831. Panel a) shows how the statistic xi varies as 
a function of the input value of a used to build the synthetic 
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CMDs. Each curve represents the run of x^ (normalised by 
its maximum value, Xj max) with a at a given cluster ring. 
The effect of mass segregation is already clear in this fig- 
ure, as the minimum value Xj min occurs at higher a as the 
outer annuli are considered, reflecting a gradual steepening 
in the PDMF. 

In order to determine a and its associated uncertainty 
from this method, 100 realisations of the model that yielded 
X7 min were run, from which a dispersion in Xy for this 
model, a-)^, was determined. Therefore, the models that sat- 
isfy the criterion 



A7 — X'y ,niin ' ^X 

are statistically of similar quality as the one that yields the 
minimum. The average value of a for these models is then 
considered as the one that best describes the PDMF, and 
the standard deviation around this average represents its 
uncertainty. Finally, panel b) in Fig. shows the a values 
determined for all annuli, again revealing variations with 
position inside the cluster. 

Following this procedure for all clusters, we built Fig. 
|S1 the CMD method counterpart of Fig. |H| Again all clus- 
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ters show evidence of mass segregation. After taking the 
uncertainties into account, these results are in very good 
agreement with those from the LP method. 



4. Discussion 

As shown in the previous section the spatial dependence 
of a clearly reveals the presence of mass segregation in all 
clusters. Except in some cases where uncorrected for incom- 
pleteness effects at lower masses may be strongly affecting 
the PDMF slope, the resufts in both LP and CMD methods 
used here are in clear agreement and within the uncertain- 
ties. 

Previous, detailed spatial determinations of the PDMF 
slope were made only for NGC 1805 (de Grijs et al. 2002a) 
and NGC 1818 (de Grijs et al. 2002a, Gouliermis et al. 
2005) . In both cases these authors apply the LP method to 
their HST/WPPC2 data to determine the P (=-a-f 1) be- 
haviour with distance to the cluster centre. In general our 
results agree with theirs, even in the case where the data 
are restricted to a higher stellar mass range (m ^ 2.0Mo) 
(Gouliermis et al. 2005). For the other clusters, there are no 
analyses of similar quality and mass range in the literature. 

Notice that a ~ 2.35 in the interval 2 < R/Rcorc < 3 
for all clusters, except for Hodge 14, for which a ^ 2.35 at 
R/Rcorc — 4. This range is close to the clusters' half-light 
(or mass) radii. 

4.1. The unresolved binarity effect 

The major systematic effect over the usual LP method to 
determine the PDMF slope is the one caused by unresolved 
binaries; a flatter PDMF (lower a) is likely to be recovered 
if binaries are not accounted for (see the discussion in Sect. 



Fig. 9. Control experiment to determine the effect of un- 
resolved binaries in the PDMF slope recovered from a sys- 
temic LP. In the vertical axis we show the difference be- 
tween the recovered slope aout and the input slope ain- 
Solid fine (dashed fine): LP generated with /bin = 0% 
(/bin = 100%). 

3). In order to assess this issue directly in our approaches, 
we generate control experiments where we simulate artifi- 
cial CMDs with /bin = 100% and an input PDMF slope 
(ain) and then recover the PDMF slope (aout) by apply- 
ing the LP method. As discussed in Sect. 3, although the 
primary and secondary stellar masses were drawn from the 
same PDMF, only primary stars with mass ^ 0.9 Mq were 
directly taken into account, in accordance with the fainter 
magnitude limit in the observed data. Furthermore, ran- 
dom pairing and a = 1.30 for 0.08 < tti/Mq < 0.8 (unseen 
mass regime) yield /bin,cff {q > 0.60) ~ 0.20. The results 
of these experiments are shown in Fig. |^ As expected, the 
recovered PDMF is flatter than the input one, and the am- 
plitude of this effect is greater for flatter input PDMFs. 
Except for ain < 1-2, however the amplitude of aout — ain 
is smaller than the uncertainties in aout- Thus, for an unre- 
solved binary fraction consistent with the one used in our 
CMD modelling process, its effect in the derived PDMF 
seems not to play a decisive role. 

4.2. Dynamical mass segregation inside the core? 

To get some clues about the nature of the mass segrega- 
tion, we evaluated the two-body relaxation time using the 
following expression (Binney & Tremaine 1987): 



trl — 



(- 



M 



6.5x10** 



ln(0.4N) ViO^Mq 



1/2 



Mr, 



3/2 



where N and M are the total number of stars and the 
total mass inside some radius R, respectively, and m, is 
a characteristic stellar mass. This is the time needed for 
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Fig. 8. The PDMF slope vs. R in the 2D CMD analysis. The lines are the same as in the Fig.|H| 



stellar encounters to redistribute the stellar kinetic en- 
ergy so that the the velocity distribution is approximately 
Maxwellian. However, it is important to keep in mind that 
this expression comes from analytical considerations. The 
tri obtained here allows a simple estimate of dynamical age 
for each cluster, since a more realistic treatment is only 
obtained through N-body simulations (Kroupa, Aarseth, 
Hurley 2001; Baumgardt & Makino 2003). 

As discussed by Stolte et al. (2002), using the present 
physical conditions of a cluster can only lead to the cur- 
rent dynamical status. We therefore decided to evaluate 
tri for approximated initial conditions, where we take the 
IMF as that from Kroupa (2002) {a = 2.30 for 0.50 < 



21.0 < V555 < 20.0 (1.5 <m/M0 

NGC 1818, NGC 1831, and NGC 1868, and 22.0 < F555 



to/Mq < 120 and a = 1.30 for 0.08 < m < 0.50Mq) 
normalised by the number of observed stars in a bright 
magnitude (massive) range. This range corresponds to 

< 2.0) for NGC 1805, 
< 

21.0 (1.2 < m/Mo < 1.5) for Hodge 14. The characteris- 
tic stellar mass to* was taken to be the median mass for 
the adopted IMF inside the observed magnitude range, 
TO, — 1.40 Mq. Table |2] shows, for all clusters, the physi- 
cal parameters estimated using this approach for all stars 
(0.08 < ui/Mq < 120) inside the core radius. The uncer- 
tainties in N (Col. 3), M (Col. 4), and t^ (Col. 5) result 
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Table 2. Dynamical parameters and timescales estimated 
for the region inside the core radius. 



Cluster 


Nobs 


iV M 
(xlO^) (xIO^Mq) 


trl 

{Myr) 


r/tri 


NGC 1805 


66 


3.9 ±0.5 2.2 ±0.3 


16 ±4 


0.64 ±0.10 


NGC1818 


140 


7.5 ±0.7 4.1 ±0.4 


48 ±7 


0.51 ±0.10 


NGC 1831 


464 


38.4 ±1.9 20.9 ±1.2 


222 ± 19 


2.26 ±0.10 


NGC 1868 


182 


16.5 ±1.3 8.9 ±0.8 


35 ±5 


25.43 ±1.00 


Hodge 14 


53 


3.7 ±0.5 2.1 ±0.3 


24 ±5 


76.25 ± 4.00 



from the statistical fluctuations in the observed star counts 
{Nobs) (Col. 2) used in the IMF normalization. 

In this table we also present the ratio t /tri in order to 
give an idea of the dynamical age of each cluster. The uncer- 
tainties in this ratio were propagated from the correspond- 
ing uncertainties in r and tri- This ratio spans a wide range 
of values (0.2 to 30). The relation between this parameter, 
evaluated for the core radius, and the corresponding PDMF 
slope (determined only by the LF method due to low statis- 
tics) is plotted in Fig. E| Although strongly based on the 
NGC 1868 result, the relation suggests that dynamically 
younger clusters tend to have steeper PDMF than dynam- 
ically older ones, as expected from the effects of dynamical 
mass segregation. We fit a linear relation (dashed line) be- 
tween these parameters (correlation coefficient (c.c.)=0.64), 
with slope dacorc/dlog(r/t,i) = — 0.38±0.30. Also note that 
all PDMF are flatter than Salpeter (marked as a horizontal 
dotted line in this flgure). By our fit one could expect that 
only clusters with \og{T/tri) ~ — 1.5 would have a PDMF 
with a Salpeter slope within their cores. 

On the other hand, if NGC 1868 is omitted from 
the linear relation (c.c.=0.47), the slope becomes flatter 
(—0.20 ± 0.27). The lack of a strong trend in the central 
a value with age could mean that mass segregation takes 
place at the onset of star formation within a cluster (pri- 
mordial segregation) or on such a short timescale that it 
leaves little room for subsequent mass segregation due to 
dynamical effects on longer timescales. There is no spe- 
cific work applying N-body simulations to address how fast 
the PMDF evolves in a typical globular cluster centre; in 
less dense environments, Kroupa, Aarseth, & Hurley (2001) 
have shown that the mass segregation in the ONC can be 
explained as a dynamical effect. Notice also that there is a 
clear PDMF flattening with dynamical age inside the core 
radius of open clusters (Bonatto & Bica 2005). 

4.3. Dynamical origin for tlie loss of lower- mass stars? 

In order to investigate a potential dynamical origin of the 
loss of lower-mass stars, we evaluated the same quantities 
as in the previous section, but considering all stars in each 
cluster. The results are presented Table|21 Although tri does 
not correspond exactly to the timescale for stellar evap- 
oration, it is expected to scale with this one (Binney & 
Tremaine 1987). As one can see in the last column of Table 
O the r/tti values again span a wide range (4x10^"^ to 1.2). 
The relation between global a (from LF and CMD meth- 
ods) and global log{T/tti) is plotted in Fig. m For both 
methods there is a trend in the sense that dynamically older 
clusters again have lower values of a. Using the results of 
the LF method, we again fit a linear relation (short-dashed 
line, c.c.= 0.87) between these parameters and find a slope 




0.5 1 

log(T/t,,) 

Fig. 10. The PDMF slope vs. log(T/tri) relation inside the 
core radius. 

Table 3. Estimates of the global dynamical parameters and 
timescales for the clusters. 

Cluster Hobs ntot rritot Ui r/iri 

(xlO'^) (xlO^Ma) (Gyr) 

NGC1805 471 30.7 ± 1.7 16.8 ± 0.9 1.1 ± 0.1 0.010 ± 0.002 

NGC 1818 774 46.3 ± 2.1 25.2 ± 1.3 1.4 ±0.1 0.018 ± 0.002 

NGC 1831 1539 109.4 ± 3.5 59.5 ±2.1 3.7 ± 0.2 0.135 ±0.010 

NGC 1868 929 63.9 ± 1.3 34.8 ± 1.5 2.0 ±0.1 0.442 ± 0.020 

Hodge 14 317 22.9 ± 1.5 12.5 ± 0.8 0.6 ± 0.1 3.073 ± 0.150 



dagiobai/dlog(r/tri) = —0.37 ±0.12. We also fit a hnear re- 
lation using the results of the CMD method (long-dashed 
line, c.c.=0.50) and find a slope of —0.13 ± 0.13. 

It is very interesting to notice that this trend is con- 
sistent with recent results in the literature for a variety of 
systems. Using N-body simulations, Baumgardt & Makino 
(2003) reproduce the observed shallowing in the global 
PDMF of globular clusters and attribute this effect to the 
lower-mass stars preferentially depleted from the cluster 
due to dynamical mass segregation. They also find that the 
details of this process are nearly independent of the start- 
ing conditions. Bonatto & Bica (2005), analysing 2MASS 
data for 11 open clusters, clearly show that these systems 
also have the same effect. In both cases the MF flattening 
reaches at least 0.8 in slope, again consistent with what 
we observe in the rich LMC clusters. Durgapal & Pandey 
(2001) show that intermediate age and old open clusters 
tend to have a smaller ratio between their present radius 
and their limiting radius as they become older. They at- 
tribute this to the effect of loss of stars. 



5. Summary and conclusions 

We have analysed the shape of the PDMF at different radii 
for a sample of 5 rich LMC clusters. The data are the re- 
sult of deep imaging with HST/WFPC2. The PDMFs were 
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Fig. 11. The PDMF slope vs. log(T/tri) for the global clus- 
ter. Open (filled) circles correspond to LF (CMD) method. 
The best linear fit for the LF (CMD) method is shown as a 
short (long)-dashed line. The dotted line corresponds to the 
Salpeter slope. 



determined by two distinct approaches: 1) a traditional 
method of converting stellar luminosities into masses; 2) a 
full modelling of the colour-magnitude diagrams of the clus- 
ters, accounting for the effects of unresolved binaries and 
photometric errors. Our results hold for both approaches 
and are insensitive to the details of how the PDMF and 
its slope are determined. Control experiments reveal that 
the unresolved binaries fiatten the PDMF recovered by the 
traditional methods. However this effect is less than the 
uncertainties in the a determination. 

We found significant mass segregation in all of them. 
The effect is expected in the sense that the PDMF is steeper 
further out than in the core: a ^ 1.80 for R < 1 Rcorc and 
a ^ Salpeter (=2.35) inside R ~ 2 ^ 3 Rcore (except for 
Hodge 14, where it occurs at i? ~ 4Rcore)- Since the global 
PDMF is also near Salpeter, we confirm previous claims 
that the PDMF evaluated around the half-mass radius is 
consistent with the global PDMF (Kroupa 2002). 

The spatial dependence of the PDMF slope was previ- 
ously presented by de Grijs et al (2002a) for the youngest 
two clusters in our sample, NGC 1805 and NGC 1818, and 
by Gouhermis et al. (2005) for NGC 1818. 

We have investigated the origin of mass segregation and 
the loss of lower-mass stars by inferring a dynamical age 
(T/tri) for each cluster. This was done in two regions: inside 
the cluster core and for the entire system. In both cases 
we notice that the dynamically older clusters (with larger 
r/tri values) tend to have shallower PDMFs. Although this 
result is based in only 5 clusters, we interpret this observed 
flattening trend as a dynamical effect, as also suggested 
by previous works for other systems (Kroupa, Aarseth & 
Tout 2001; Baumgardt & Makino 2003; Bonatto & Bica 
2005). In fact, it is hard not to expect some dynamical 



signature in the PDMFs for those LMC clusters with quite 
different dynamical ages. On the other hand, the PDMF 
flattening (that can reach Aa ^ 1.0) could be interpreted 
as primordial mass segregation, at least in the case of the 
central regions. Or perhaps dynamical mass segregation in 
the LMC cluster's cores occurs fast enough that any age 
dependence of the PDMF would be seen only in a sample 
of young clusters. Again we stress the importance of N- 
body simulations that configure a unique tool for providing 
reliable answers. 

Although a scatter of ^^ 1.0 in a is expected due to 
Poisson noise and the dynamical evolution of stellar clus- 
ters (Kroupa 2001), the aforementioned trends are in accor- 
dance with what one expects from dynamical arguments. In 
order to improve this dynamical investigation, and there- 
fore confirm or reject these trends, we are applying our 
PDMF determination techniques to a larger number of 
CMDs of LMC/SMC clusters (Brocato et al. 2001). 
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